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Fission of doubly charged metal clusters is studied using the open-shell two-center 
deformed jellium Hartree-Fock model and Local Density Approximation. Results of 
calculations of the electronic structure and fission barriers for the symmetric and 
asymmetric channels associated with the following processes 

Na^g — > Naf 5 + Na^ and Na^ — > 2Na§ are presented. The role of the exact 
exchange and many-body correlation effects in metal clusters fission is analysed. It 
is demonstrated that the influence of many-electron correlation effects on the height 
of the fission barrier is more profound if the barrier arises nearby or beyond the 
scission point. The importance of cluster deformation effects in the fission process is 
elucidated with the use of the overlapping-spheroids shape parametrization allowing 
one an independent variation of deformations in the parent and daughter clusters. 



I. INTRODUCTION 

Fission of charged atomic clusters occurs when repulsive Coulomb forces, arising due to 
the excessive charge, overcome the electronic binding energy of the cluster jl], 0, H]. This 
mechanism of the cluster fission is in a great deal similar to the nuclear fission phenomena. 
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Experimentally, multiply charged metal clusters can be observed in the mass spectra when 
their size exceeds the critical size of stability, which depends on the metal species and cluster 
charge [|, j|, || [7|, ||, [|. For clusters above critical size, simple evaporation of neutral species 
is the dominant fragmentation channel, while below the critical size, fission into two charged 
fragments is more probable. At a low temperature fissioning clusters can be metastable 
above a certain size, because of the existence of a fission barrier. 

Initially, theoretical studies of cluster stability were based on pure energetic criteria that 
only involved the energies of the initial and the final states [10, 11, 12|| . Later, a simple 
one center Liquid Drop Model (LDM), which was initially suggested by Lord Rayleigh in 



1882 JT3| and later widely used in nuclear physics, was adapted to charged metal clusters 



|LH1 . In this model, one introduces the "fissility parameter", X = E^^f/2E s s p ^ r e p, which is 
proportional to the ratio of the Coulomb to surface energy of charged spherical liquid drop 
|T3| . The fissility parameter distinguishes the situations when cluster is unstable, metastable, 
or stable. The investigation of the Rayleigh instabilities in multiply charged sodium clusters 



has been done in |15|], where reasonable agreement with experimental data was found. 

In spite of the fact that the simple LDM qualitatively describes the fission process it fails 
to reproduce experimental data in full detail. This happens because the LDM does not take 
into account shell effects. It has been shown that the shell effects are important in nuclear 
fission ]TE[ and even more important in fission of metal clusters [0. One can describe the 
shell effects in metal clusters using the Shell Correction Method (SCM), originally developed 



in nuclear physics |16|, [17], |18| . This method was adapted for metal clusters in [|19], ^(J, [21], 



H H HI- 



The asymmetric two-center-oscillator shell model (ATCOSM), introduced in [25| for nu- 
clear fission is quite successful in prediction of the fission barriers. This model was also 
applied for the description of metal cluster fission [2(|. Although ATCOSM method uses 
single electron model potential it has the significant advantage in comparison to other mod- 
els allowing one simple shape parametrization and an independent variation of deformations 



in the parent and daughter clusters [26 



The microscopic description of energetics and dynamics of metal cluster fission process 



based on molecular dynamic (MD) simulations has been performed in |57], |28], ^9] [3(]] using 
the local-spin-density-functional method. This method is, however, strongly restricted by 
the cluster size, because of computational difficulties, and thus is usually applied to the 
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small metal clusters with the number of atoms N < 20. 

Fission process of metal clusters can be also simulated on the basis of the jellium model, 
which does not take into account the detailed ionic structure of the cluster core. Jellium 
model considers the electrons in the usual quantum mechanical way, but approximates the 
cluster core potential by the potential of the homogeneous positively charged background 
and, therefore, is better applicable for lager cluster sizes (see e.g. |H| for review). Most of 
the electronic structure calculations of the jellium metal clusters have been performed using 
self-consistent Kohn-Sham Local Density Approximation (LDA) f32j. The LDA jellium 
calculations for metal cluster fission can be grouped |26f into two categories according to the 
fragment shape parametrization, namely, the two-intersected-spheres jellium [[33, [34, |35|] and 



variable-necking-in parametrization |36|, [37], |38]. It has been shown that the cluster shape 



parametrization must be flexible enough to account for the majority of effects generated by 
the shell structure of the parent and daughter clusters, which in general have not spherical 



but deformed shapes [26 



The important feature of the LDA method consists in the fact that it takes into ac- 
count many-electron correlations via the phenomenological exchange-correlation potential 
(see e.g. [^, fJ(J for review). However, so far, there has not been found the unique poten- 
tial, universally applicable for different systems and conditions. As a result there is a "zoo 
of potentials" |?T| valid for special cases. These potentials, of course, do exist in principle 
as unique quantities but are not actually understood, so they cannot serve as a satisfactory 
basis for achieving a physical interpretation. 

Alternatively, one can develop direct ah initio methods for the description of electronic 
properties of metal clusters. It can be achieved by using the Hartree-Fock (HF) approxi- 
mation and by the construction on its basis the systematic many-body theories such as the 



random phase approximation with exchange [42|, many-body perturbation theory or the 



Dyson equation method [43], [44], [45]] . Based on fundamental physical principles these models 
can be refined by extending the quality of the approximations, while the physical meaning 
of the effects included are clearly demonstrated and thus give more accurate characteristics 
of metal clusters than LDA. 

Originally, the Hartree-Fock model for the metal cluster electron structure has been 
worked out in the framework of spherically symmetric jellium approximation in [JB, H7J. It 
is valid for the clusters with closed electronic shells that correspond to magic numbers (8, 
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20, 34, 40,...). On the basis of the Hartree-Fock approach the dynamic jellium model has 
been proposed |^8|, [0| . This model treats simultaneously the vibrational modes of the ionic 



jellium background, the quantized electron motion and the interaction between the electronic 
and the ionic subsystems. In particular, the dynamic jellium model allows to describe 
the widths of electron excitations in metal clusters beyond the adiabatic approximation. 
The open-shell two-center jellium Hartree-Fock approximation valid for metal clusters with 
arbitrary number of the valence electrons has been developed in J5(| [51]]. The two-center 



jellium HF method treats the quantized electron motion in the field of the spheroidal ionic 
jellium background in the spheroidal coordinates. This method has been generalized and 



adopted to study of the metal clusters fission process in our recent work |p2| , where barrier 
for the symmetric fission cannel Na%g — > 2Na^ was calculated. 

In the present work we investigate the role of the exchange and correlation effects in 
metal cluster fission process on the basis of both the Hartree-Fock and LDA methods. 
Both symmetric and asymmetric fission channels for the Na\^ and Naf^ parent clusters 
are considered. Comparison of results of the two approaches allows us to illustrate the 
importance of the exchange component of the many-electron interaction in the fission process 
and make important conclusion about the relative role of the two different channels of the 
reaction. 

The atomic system of units, \e\ — m e — h — 1, has been used throughout the paper, 
unless other units not-indicated. 

II. THEORETICAL METHODS 

A. Two-center jellium model and cluster shape parametrization 

According to the main postulate of the jellium model, the electron motion in a metal- 
lic cluster takes place in the field of the uniform positive charge distribution of the ionic 
background. For the parametrization of the ionic background during the fission process we 
consider the model in which the initial parent cluster having the form of the ellipsoid of 



revolution (spheroid) splits into two independently deformed spheroids of smaller size [52 



The two principal diameters ak and bk of the spheroids can be expressed via the deformation 
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2 + 8 k \ 2/ \ t /2-4 xV3 



parameter 8 k as 

ak= {j^f k ) Rk > bk = [jTT k ) Rk - (1) 

Here partial indexes k = 0, 1, 2 correspond to the parent cluster (A; = 0) and the two daughter 
fragments (k = 1,2), i?fc (fc = 0, 1, 2) are the radii of the corresponding undeformed spherical 
cluster. The deformation parameters 8 k characterize the families of the prolate (8 k > 0) and 
the oblate (8 k < 0) spheroids of equal volume V k = ^Ka k b\j3 = 4nRl/3. 

The radii of the parent and the resulting non overlapping daughter fragments are equal to 

1 /3 

R k = r s N k , where N k is the number of atoms in the fcth cluster, and r s is the Wigner-Seitz 
radius. For sodium clusters, r s = 4.0, which corresponds to the density of the bulk sodium. 
For overlapping region the radii R\(d) and R2(d) are functions of the distance d between the 
centers of mass of the two fragments. They are so determined that the total volume inside 
the two spheroids equals the volume of the parent cluster 47tRq/3. 

The ions charge density p(r) is kept uniform including the overlapping-spheroids region, 

p c , ( x i + yi)/bl + (z + d/2) 2 /al<l 

p c , {x 2 + y 2 )/bl + {z-d/2f/al<l (2) 
0, otherwise. 

Here p c = Zq/Vq is the ionic charge density inside the cluster and Z is the total charge of 
the ionic core. 

The electrostatic potential U core (r) of the ionic background can be determined from the 
solution of the corresponding Poisson equation: 

AU core (r) = -47rp(r). (3) 



p(r) 



B. Hartree-Fock and LDA formalism 

The Hartree-Fock equations can be written out explicitly in the form (see, e.g., H53|): 

(- A/2 + U core + U H f) I a > = e a | a > . (4) 

The first term here represents the kinetic energy of electron a, and U core its attraction to 
the cluster core. The Hartree-Fock potential Uhf represents the average Coulomb interac- 
tion of electron a with the other electrons in the cluster, including the non-local exchange 
interaction, and e a describes the single electron energy. 
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According to the Density Functional Theory (DFT), the ground state energy is a mini- 



mum for the exact density of the functional of the density of the system . A self-consistent 
method for calculation of the electronic states of many-electron systems was proposed by 
Kohn and Sham This method leads to the Kohn-Sham LDA self consistent equations, 



which are similar to the Hartree equations: 

(-A/2 + U core + U H + V xc ) | a > = e a | a > . (5) 

Here Uh is the Hartree potential, which represents the direct Coulomb interaction of elec- 
tron a with the other electrons in the cluster, but does not take into account the non- 
local exchange effects, while V xc is the phenomenological density dependent local exchange- 



correlation potential. In the present work we use the Gunnarsson and Lundqvist model |55 
for the LDA electron exchange-correlation energy density e xc , which reads as 

e xc (Pei(r)) = -\ (J^j ' ^ " °- 0333 G (r.W/H.4) • (6) 

Here r s (r) = (3/47rp e ;(r)) 1//3 is a local Wigner-Seitz radius, while p e ;( r ) i s the electron density 
in the cluster, and the function G(x) is defined by following relation: 

G(x) = (l + a; 3 )ln(l + ^-x 2 + |-i (7) 

The exchange-correlation energy density e xc (p e i(r)), defines the LDA exchange-correlation 
potential V xc (p e i{r)) as 

= *Mr)UMr))] = 
opei{r) 

9 \ 1/3 1 „ „„„„ , A 11.4 



0.0333 In 1 + 



4tt 2 J r s (r) \ r s (r 

The Hartree-Fock (Hj) and LDA (|^) equations have been solved in the system of the 
prolate spheroidal coordinates as a system of coupled two-dimensional second order partial 
differential equations. The partial differential equations have been discretized on a two- 
dimensional grid and the resulting system of linear equations has been solved numerically 



by the successive overrelaxation method [56]. This technique is different from that we 



have used in our previous works |50| , |5T| , |52fl , where the expansion of wave functions and 
potentials over spheroidal harmonics in the prolate and oblate spheroidal coordinates in 
one dimension and a numerical expansion for the second dimension have been carried out. 
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The third dimension, the azimuthal angle, has been treated analytically in both methods. 
Our calculations show that the partial-expansion method is effective for slightly deformed 
systems, for which only few terms in the expansion are necessary to take into account in 
order to achieve sufficient accuracy. However, for strongly deformed systems, or process like 
fission, the direct two-dimensional integration is more efficient. 

The important characteristic of the cluster, which defines its stability is the total energy 
E tot . The total energy of the cluster is equal to the sum of the electrostatic energy of the 
ionic core E core and the energy of the valence electrons E e i. 

E t ot = E core + E d . (9) 

The electrostatic energy of the cluster ionic core is equal to 

E CO re = ~ J p(r)[/ core (r)dr. (10) 

In the Hartree-Fock approximation, the electronic energy E e \ is given by the general 
expression [Q: 

E el F = J2< a \ ~ A /2 + U core \a> + 
a 

- qaqb [c( abk )E k {a, b) + d{abk)G k {a, b)} , (11) 

abk 

where a and b run over all shells. The values F k (a,b) and G k (a,b) in the Eq. (01]) are the 
Coulomb and exchange Slater integrals, q a and are the occupation numbers for orbitals 
a and b, respectively. The Hatree-Fock coefficients c(abk) and d(abk) for the Coulomb and 
exchange energy contributions depend on the occupation numbers (see for details |53|). 



In the framework of LDA the electronic energy of the system is given by |54 

E 



E el DA = V < a I "A/ 2 + Ucore \ > + 



1 J P^) ^ + J p^eUPe^dr, (12) 
where the latter term represents the exchange-correlation energy. 



III. NUMERICAL RESULTS 



Let us present and discuss the results of calculations performed in the model described 
above. We start our consideration with the analysis of electronic configurations and occupa- 



8 




FIG. 1: Hartree-Fock single electron energy levels for the asymmetric fission channel Na^ — > 
Naj + Na$ as a function of fragments separation distance d, when both the parent and daughter 
clusters are spherical (So = 5% = 82 = 0). The evolution of the cluster shape during the fission 
process is shown on top of the figure. Horizontal lines on the right hand side of the figure mark 
the Hartree-Fock energy levels for the free daughter fragments. 

tion numbers alterations during the fission processes of the doubly charged sodium clusters 
Na^Q and Na^, which we perform using the two overlapping sphere parametrization model. 
The second part of our discussion is devoted to the cluster energetics and formation of fis- 
sion barriers, for the fission channels considered. We compare the results obtained for the 
two overlapping spheres and the two overlapping spheroids parametrization, as well as the 
variable-necking-in type of the shape parametrization and demonstrate the crucial impor- 
tance of the cluster shape deformations in the fission process. 

Figures [l|-|3| show the Hartree-Fock single-electron energies Si as a function of the fragment 
separation distance d for the following asymmetric Na\$ — > Ndj + Na^, Nal$ — > Naf 5 + 
Na^ and symmetric Na^g — > 2Na^ channels respectively, when both parent and daughter 
clusters are assumed to be spherical, and, hence, for any d the deformation parameters are 
equal to zero 8 = 81 = <5 2 = 0. The snapshots of the cluster shape evolution during the 
fission are shown on top of each figure. 

The initial electronic configuration for the parent Na^ cluster is lcr 2 2cr 2 l7r 4 , where 2al% 
level is degenerated, due to spherical symmetry of the system with closed electronic shells. 
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FIG. 2: The same as figure |l] but for asymmetric fission channel Nalg — > Naf 5 + Na^ . Lowest 
unoccupied states are shown by dashed (Act) and short dashed (15) lines. 
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FIG. 3: The same as figure |l| but for symmetric fission channel Na^ — ► 2Na^ . Lowest unoccupied 
states are shown by dashed (4a) and short dashed (2ir) lines. 

With increasing fragments separation distance d, the spherical symmetry of the parent 
cluster breaks down, that results in splitting of the single electron energy levels according to 
projection of the angular momentum along the z-axis. A similar splitting of the energy levels 
has been studied for deformed nuclei (see e.g. |57| ) and clusters f5(J, [5l|]. For high enough 
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separation distances beyond the scission point, indicated by vertical arrow A in figure [I], 
energy levels approach to their limiting values marked by horizontal lines on the right hand 
side of the figure, being the electron energy levels of the free fission fragments Na^ (la 2 l7r 4 ) 
and Naj (la 2 ). 

The electronic configuration and the orbital occupation numbers for iVa 2 ^ cluster exhibits 
several alterations during the fission process, as it is shown in figures |2] and [3[ This happens 
due to the fact that different electronic configurations minimize the total energy of the 
cluster at different separation distances. The parent cluster iVa 2 ^" has the initial electronic 
configuration lcr 2 2cr 2 l7r 4 27r 2 3cr 2 15 4 . Following to the asymmetric channel of fission Na1$ — > 
Naf 5 + Na^ (see figure |]) the new configuration la 2 2a 2 l7r 4 3a 2 4a 2 l<5 4 becomes preferable, 
when two electrons transfer from the half-filled 2tt state to initially unoccupied 4a state. 
This happens for the separation distance d > 14.2 a.u. (marked in figure |^ by the solid 
vertical arrow A) before the scission point (vertical arrow B) d = 15.3 a.u. The order of 
the energy levels manifests several re-distributions during the fission process, and finally for 
high enough separation distance (d > 25 a.u.), its energy levels correspond to the free fission 
fragments Naf 5 (la 2 2a 2 l7T 4 3a 2 15 4 ) and Na+ (la 2 ). 

For the symmetric fission channel iVa 2 ^ — > 2Na^ (figure |3|) the energy levels show even 
more complicated behaviour. Thus, for the separation distance d > 0.7 a.u. (marked in 
figure [5] by solid vertical arrow A) the new intermediate configuration Ia 2 2a 2 l7r 4 l5 2 27r 4 3a 2 
becomes preferable, when two electrons transfer from the occupied 15 to the half-filled 2ir 
state. Further increasing of the separation distance, leads to the transition of the two residual 
15 electrons to the initially unoccupied 4a state for d > 9.25 a.u. (marked in figure |3| by 
the solid vertical arrow B) forming the final electronic configuration la 2 2a 2 l7r 4 3a 2 2vr 4 4a 2 
before the scission point (vertical arrow C) d — 16.64 a.u. The order of the energy levels 
manifests several re-distributions during the fission process, and finally for high enough 
separation distance (d > 25 a.u.), its ordering corresponds to la2aln which determined by 
the magic spherical Nag products. They are doubly degenerated as compared to the initial 
configuration, since there are two Na^ fragments in the system. Such a behaviour of the HF 
levels as a function of d is quite similar to those following from the ATCOSM simulations 

Figure ^ presents fission barriers for the asymmetric channel Na^ —>■ Na^ + Na^ as a 
function of the fragments separation distance d. In order to perform the accurate comparison 
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FIG. 4: Fission barriers in the two-center deformed jellium Hartree-Fock (solid lines) and LDA 
(dashed lines) approaches calculated in this work for the asymmetric channel No^q — > Na^ + Na^ 
as a function of fragments separation distance d. In (a), both parent and daughter clusters are 
spherical, 5q = 5\ = 62 = 0. In (b), deformations of the parent and daughter clusters are taken 
into account. The zero of energy put at d = 0. The evolution of cluster shape during the fission 
process is shown on top of figures (a) and (b) for both models. We compare our results with those 
derived in ATCOSM (dash-dotted line) f26j. 
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of fission barriers derived in the two-center deformed jellium HF and LDA models with the 



ATCOSM results Pq| , we have used two different type of shape parametrization. Thus, 
the upper part in figure |] shows the barriers for fission of a spherical parent cluster into 
two spherical daughter fragments (i.e. 5q = S\ = 62 = in our model). This type of 
parametrization, known as the two-intersected spheres parametrization, was used in many 
works |33| , P^ , [35]. The low part in figure ^ shows fission barriers derived on the basis of 



parametrization accounting for an independent deformation of parent and daughter clusters 
in order to minimize the total energy of the system for any distance d. The evolution of 
cluster shape during the fission process is shown for the HF and LDA models on tops of the 
corresponding figures. Note, that the variable-necking-in parametrization has been used in 



the ATCOSM calculation [26 



Solid lines in figure |] are the result of the two-center jellium HF model, while dashed 
curves have been calculated in LDA. Dash-dotted lines show the ATCOSM barriers calcu- 
lated in |p6| . The zero of energy put at d = 0. 

Figure [| (a) demonstrates a good agreement of the HF and ATCOSM fission barriers 
heights Ahf = Aatcosm = 0.95 eV. The LDA value for the fission barrier, Alda = 1-30 
eV, slightly exceeds the HF and ATCOSM ones. In HF and LDA schemes the fission barrier 
maximum is located just behind the scission point (marked by vertical arrow A). 

The two-intersected spheres parametrization is not adequate for the description of the 
process Na^ — > Na^ + Na^, because Na^ daughter fragment has an open-shell electronic 
structure, and therefore its shape is not spherical because it undergoes Jahn- Teller distortion 



24j . The deformed jellium Hartree-Fock and LDA calculations show that the ground state of 



the Na^ cluster has the oblate shape with the deformation parameter 8\ = —0.68. This value 
is in a good agreement with other theoretical estimates ||50|| . The oblate shape deformation 
of the daughter cluster reduces the final total energy of the system E tot by —1.32 eV for 
the HF, and by —1.05 eV for the LDA model calculations, in comparison with the spherical 
case. 

Part (b) , figure |] shows fission barriers for the asymmetric channel Na^ — > Na^ + Na^ , 
when spheroidal deformations of the parent and daughter clusters are taken into account. 
We have minimized the total energy of the system over the parent and daughter fragments 
deformations with the aim to find the fission pathway corresponding to the minimum of the 
fission barrier. We have also used the assumption of continuous shape deformation during 
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the fission process. The deformation of the cluster fragments changes drastically the fission 
energetics in comparison with what follows from the two-intersected spheres model. In 
the framework of the two-center deformed jellium Hartree-Fock approximation, the parent 
cluster NclIq becomes unstable towards the asymmetric channel No^q — > Na^ + Na^. 

The LDA simulations with deformation of the daughter fission fragments take into account 
the decrease of the total energy of the system (see figure |]b). In particular, this results in the 
appearance of the local minimum in the energy curve at d — 7.2 a.u., corresponding to the 
formation of the super deformed asymmetric prolate shape of the parent Naf^ cluster before 
the scission point A. The latter is located at d — 10.4 a.u. The allowance for deformation 
of the parent cluster and the fission fragments reduces the LDA fission barrier up to the 
value Alda = 0.16 eV, which is in rather poor agreement with the result of ATCOSM 
Aatcosm = 0.52 eV. 

It is necessary to note that the shape of the cluster fragments in ATCOSM has been 



parametrized by two spheroids of revolution connected by a smooth neck p6| . Our calcu- 
lations show that the oblate shape of Na^ fragment is formed at the initial stage of fission 
process, for separation distances before the scission point. Moreover, in the vicinity of the 
scission point, where the interaction between the two daughter fragments Na^ and Na% 
is strong, the oblate Na^ fragment is even more deformed than a free one. This means 
that it is more favourable for two fragments to split at shorter distances rather than to 
be connected by a smooth neck, making the system more prolate. This fact explains why 
using of necking type of shape parametrization, leads to the higher barrier as compare to 
our parametrization. 

Comparison of the asymmetric Na\% — > Naf 5 + Na^ and symmetric Na1g — > 2Na,Q 
fission channels of the parent Naf^ cluster is a subject of particular interest, because there 
have to be a competition between these two channels involving magic cluster ions Na^ 
and Nag. In || 23), |27] it was noticed that namely in fission of the Na±g cluster a magic 
fragment other than Na^ becomes the favoured channel. Figure [| shows fission barriers for 
the asymmetric channel Na%g — > Naf 5 +Na^. We have started from the initial configuration 
corresponding to oblate shape of the parent Naf^ cluster with the deformation parameter 
5 = —0.35. The oblate deformation reduces the total energy of the cluster A"a^g by —0.58 
eV for the HF and by —0.48 eV for the LDA model in comparison with the spherical case. 
We have minimized the total energy of the system over the deformation parameters of the 
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FIG. 5: The same as figure but for the asymmetric channel Na^ — > Naf^ + Na^ 

parent and daughter fragments during the fission process for any separation distance d. The 
evolution of the fragment shapes is shown on top of the figure for both HF and LDA models. 
The daughter fragment Naf 5 has an oblate shape with deformation parameter 8\ = —0.6, 
while Na% is spherical, i.e. 82 = 0. 

The total energy as a function of the fragment separation distance has a maximum 
(marked by vertical arrow A for HF, and A for LDA), arising due to the alteration of 
the electronic configuration Icr 2 2cr 2 l7r 4 27r 2 3er 2 lr5 4 — > l(j 2 2<7 2 l7r 4 3<7 2 4cr 2 l<5 4 . These maxima 
on the energy curves define the fission barrier hights, being equal to Ahf = 0.36 eV for 
Hartree-Fock and &lda = 0.50 eV for LDA. It is interesting to notice that the LDA total 
energy curve has a pronounced minimum at d — 12.5 a.u., located beyond the scission point 
d = 11.1 a.u. We have marked the scission point by vertical arrow B, both for HF and 
LDA. This minimum means that a quasistable state of the supermolecule Naf 5 + Na% can 
be created during the fission process. However, the appearance of the minimum and thus 
the stability of the super molecule is rather sensitive to the model chosen for the description 
of exchange and correlation inter-electron interaction. This is already clear from the fact 
that such a minimum does not appear in the HF simulations (see figure ||). 

The ATCOSM model calculation gives the value of the fission barrier Aatcosm = 1-35 
eV, which is inconsistent with the HF and LDA results derived in our model. Such a differ- 
ence can be explained as a result of variable-necking-in type of the shape parametrization, 
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FIG. 6: The same as figure but for the symmetric channel Na1$ — > 2Na 



which has been used in [26]. In the case of asymmetric Na^ — > Naf 5 + Na^ channel, 
the parent as well as one of the daughter fragments have an oblate shape, therefore shape 
parametrization model with prolate-like additional neck is not natural, and results in in- 
creasing the fission barrier. 

Figure |^ shows the dependence of total energy E tot on separation distance d for the sym- 
metric channel Naf^ — ► 2Na^. The parent cluster changes its shape from oblate to prolate 
one on the initial stage of the fission process (d « 1 a.u.). This transition is accompanied 
by the first re-arrangement of the electronic configuration (marked by vertical arrow A for 
HF and A' for LDA). This process has the barrier Ahf — 0.63 eV and Alda = 0.48 eV. 
On the next stage of the reaction the prolate deformation develops resulting in the highly 
deformed cluster shape, as it is shown on top of figure |^. At the distance d ~ 11 (marked 
by vertical arrow B for HF, and B' for LDA) the electronic configuration reaches its final 
form being the same as in the spherical Nag products. In this case, the variable-necking-in 
type of fragments shape parametrization used in ATCOSM does not brake the symmetry 
of the fragments and therefore the agreement between the two overlapped spheroids HF or 
LDA models and ATCOSM variable-necking-in approach is much better than in the case 
of asymmetric fission channel. Indeed, the total fission barrier for the symmetric channel 
Nal$ — > 2Na,g is equal to Ahf = 0.63 eV and Alda = 0.48 eV in the two-center jellium 
Hartree-Fock and LDA models respectively. These values are in a good agreement with the 
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TABLE I: Summary of the fission barrier heights in (eV) calculated in this work and their com- 
parison with the results of other approaches. 



Channel 




HF 


LDA 


ATCOSM H 


MD |29| 


MD H 


MD |27| 




Na$ + Na£ 





0.16 


0.52 


0.67 


0.50 


0.71 




Na^ 5 + Na% 


0.36 


0.50 


1.35 


0.50 






Na 2 +^ 


2Na£ 


0.63 


0.48 


0.52 


0.52 







ATCOSM result A ATCO sm = 0.52 eV M 



In table|| we have summarized the results of the HF and LDA barrier heights calculations 



for the considered fission channels and compared them with the results of ATCOSM [26 



and MD simulations 27, 29, 30 



The height of the fission barrier for Na±Q cluster in the two-center deformed jellium LDA 



model is 0.51 eV lower than its value following from the MD simulations [29|. Molecular 



dynamics simulations performed in |29[ also were based on the density-functional theory, 
but included full ionic structure of the cluster. Since both methods apply the same form of 



the density functional fl55 |, the discrepancy in the fission barrier heights can be attributed 
to the manifestation of the influence of the cluster ionic structure in the fission process. One 
can expect that the influence of the detailed ionic structure has to decrease with the growth 
cluster size, making the jellium model approach more and more accurate. However, we also 
want to notice here that different schemes of MD simulations [27, ^9| lead to somewhat 
different fission barrier heights (see Table [Hl|) . 

For the Na^ cluster, we report a very good agreement of the heights of fission barriers 



derived in the jellium LDA model and MD |[29|| . It is interesting to note that MD simulations 
predict the asymmetric channel Na1$ — > Naf 5 + Na^ to be a little more favourable, however 
the jellium LDA model says in favour of the symmetric path Naff — > 2Nag. Although, in 
both models the difference in the heights of symmetric and asymmetric fission barriers is 
rather small, being equal to 0.02 eV. 

Finally, let us compare the results obtained in the HF and LDA models. The Hartree- 
Fock model takes into account the exchange inter-electron interaction the most exactly, 
but it neglects many-electron correlations, playing quite a significant role in metal cluster 
energetics. Thus, LDA calculations with the Gunnarsson-Lundqvist's exchange-correlation 
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FIG. 7: Difference between the LDA and HF total energies as a function of fragments separation 
distance for the fission channel Naf^ — > Na^ + Na^ . The evolution of the cluster shape during 
the fission process is shown on top of the figure. Vertical arrow A corresponds to the scission point. 

potential |55j lower the total energy E tot of the parent Naf^ and Naff clusters on —5.58 
eV and —12.03 eV respectively as compared to the energies calculated in the Hartree-Fock 
approximation. Such a difference between the total energies of the HF and LDA jellium 
clusters is a result of the manifestation of many-electron correlation interaction. The many- 
electron correlation interaction reduces at distances d beyond the scission point, when the 
parent cluster splits on the two fragments, due to decrease of the number of interacting 
electrons in the system. 

Figure [7| shows the difference between the LDA and HF total cluster energies E^ f A — E^f 
as a function of fragments separation distance for the fission channel No^q — > Na^ + Na^. 
For illustrative purposes we use here the two overlapping sphere parametrization to dismiss 
the effect of the clusters fragment shape alteration along the fission pathway, which can be 
different in the HF or LDA schemes of calculation. 

It is seen from figure [7] that the difference between the LDA and HF total energies 
^to? A ~ EtoF does not change significantly for separation distances below the scission point, 
in spite of the strong deformation of the parent cluster. In the narrow region of d, nearby 
the scission point the value Ef£ A - EfJ increases and it becomes constant again above 
the scission point. Such a behavior is in a great deal similar, at least qualitatively, to the 
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all considered fission channels. It remains valid even if deformations of the parent and 
daughter fragments are taken into account. This fact has a simple physical explanation. 
The many-electron correlations reduce the total energy. During the fission process the 
many-electron correlation interaction between electrons belonging the two different cluster 
fragments vanishes, which results in the increase of the total energy of the system. 

In spite of the fact that many-electron correlations reduce significantly the total energy 
of the Na\^ and Na\% jellium clusters in comparison with the HF values, the difference in 
heights of the fission barriers obtained in the HF and LDA models is only about 0.15 eV. 
From the analysis carried out above one can conclude that many-electron correlations do 
not influence significantly on the height of the fission barrier if the latter arises well below 
before the scission point. In the cases when the barrier is created nearby or above the 
scission point, accounting for many-electron correlations becomes essential. For example, 
many-electron correlations play the crucial role for the No\q — > Na^ + Na^ fission channel. 
In this case the HF model predicts qualitatively wrong barrierless scenario of the fission. 
In the contrary, the agreement between the HF and LDA models is much better for both 
symmetric and asymmetric fission channels of the Nafl cluster. Note also that the HF 
model, predicts the asymmetric fission channel Naf£ — > Naf 5 + Na^ to be more favourable 
in comparison with the symmetric one. 



IV. SUMMARY 



We have developed the open-shell two-center deformed jellium Hartree-Fock and LDA 
method for the description of metal clusters fission process. The proposed two overlapping 
sheroids shape parametrization allows one to consider independently a wide variety of shape 
deformations of parent and daughter clusters, and to investigate the role of deformation 
effects on the cluster fission process. The proposed type of shape parametrization is prefer- 
able as compared to variable-necking-in one, in particular when the fission fragments have 
an oblate shape. The role of many-electron correlation effects in metal clusters fission is 
analysed. The described Hartree-Fock model forms the basis for further systematic devel- 
opment of the more advanced ab initio many-body theories for the process of metal clusters 
fission. 
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